Domain Growth in Chiral Phase Transitions: 

Inertial Dynamics 



Awaneesh Singh 1 , Sanjay Puri 1 and Hiranmaya Mishra 1 ' 2 
1 School of Physical Sciences, Jawaharlal Nehru University, New 

Delhi-110067, India. 
2 Theory Division, Physical Research Laboratory, Navrangpura, 
Ahmedabad-380009, India. 

Abstract 

We investigate the kinetics of phase transition for chiral symme- 
try breaking after a collision with high-energy heavy ions. We use a 
O I Langevin description for domain growth kinetics in the chiral phase 

^ ■ transition. The Langevin equation of motion for the order parame- 

q_) _ ter field includes dissipation and the inertial term with a Gaussian 

random force. We study the far-from-equilibrium ordering dynamics 
subsequent to a quench from the massless quark phase to the mas- 
sive quark phase and discuss the effect of dissipation on the dynamic 
magnification. After a discussion on the growth of the bubble, we 
introduce the assumption of scale and to show how a dynamic corre- 
lation length emerges. We find a simple length scale corresponding to 
the Allen-Chan (AC) growth law i.e., L(t) ~ t 1 / 2 . 
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1 Introduction 



Strongly interacting hadronic matter is expected to undergo a phase transi- 
tion at sufficiently high temperature and/or baryon density. For vanishing 
baryon density, this is a prediction from abinitio calculations like the lattice 
simulation of the quantum chromo dynamics (QCD) while for finite baryon 
density, this is expected from different effective models of the strongly in- 
teracting system. Heavy ion collisions provide an opportunity to scruti- 
nize the phase structure of QCD experimentally. The high temperature and 
low baryon density phase was explored in the relativistic heavy ion collision 
(RHIC) experiments at Brookheaven and it will be further studied in the 
large hadron colliders (LHC) experiments. The planned experiments like 
beam energy scan (BES) at RHIC pQ, compressed baryonic matter (CBM) at 
GSI [2] and nuclotron-based ion collider facility (NICA) at Dubna [3] intend 
to study different aspects of the thermodynamic properties of QCD at finite 
chemical potential, such as the expected critical point of QCD, first order 
phase transition as well as the mixed phase structure. 

However, it is challenging to extract the thermodynamic properties of 
quark hadron phase transition from the nuclear collision experiments. Apart 
from the system size being small, the dynamics here is rather fast and can 
prevent the thermal equilibrium being established globally. Therefore, it is 
likely that the nonequilibrium effects play an important role in the evolution 
of the fireball. It may be emphasized here that an understanding of the 
equilibrium phase diagram is not sufficient to surmise the properties of such 
a system. One also has to understand the kinetic process that drive the phase 
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transition and the properties of the nonequilibrium structures that the system 
goes through on its way to reach equilibrium j3]. In this regard, both the 
critical dynamics and the far- from- equilibrium kinetics of chiral transition has 
attracted recent attention. In the study of critical dynamics (i.e., the time- 
dependent behavior in the vicinity of the critical point), much interest has 
been on the signatures of the critical end point (CEP) of QCD [5], El El IH] - In 
this paper, we focus on the second aspect: the far-from-equilibrium kinetics 
i.e., the evolution of the system after a quench from a disordered phase at 
high temperature to an ordered phase of low temperature with nonvanishing 
quark antiquark condensates. 

In this context, the relaxation to equilibrium using Langevin equation has 
been attempted in Ref . [9] , where the authors studied the early time dynam- 
ics of the spinodal decomposition and the effect of dissipation on the spinodal 
instability. Further, the bubble nucleation kinetics in chiral transition was 
studied in Ref. [10]. A time- dependent Ginzburg Landau equation (TDGL) 
was derived in Ref. [11] starting from a non-ideal, nonrelativistic hydrody- 
namics for coupled order parameters. The authors here clarified the effect 
of viscosity in the ordering kinetics. Further, in Ref. [12], the amplification 
of the spinodal fluctuation was studied within a fluid dynamical model for 
the nuclear collisions. Here, the study was mostly focused on the evolution 
in the linearized regime which showed an exponential growth of the initial 
fluctuations. 

We have recently attempted the far-from-equilibrium kinetics of chiral 
phase transition using a TDGL equation [13], which was complementary 
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to the studies of Refs. [HI QUI [HI E2] and investigated the late stages of 
phase separation kinetics in quark matter. The Ginzburg Landau free energy 
functional was obtained from an expansion of the thermodynamic potential 
of a two flavor Nambu-Jona-Lasinio (NJL) model. The nonlinearities in the 
evolution equation made the exponential growth of the initial fluctuations 
stabilize. The quenching scenario was studied where the quenching is done 
through the second order as well as the first order line in the phase diagram. 
The domain growth in both types of quenching was studied and different 
quantitative features of the coarsening morphologies were also discussed. 

Let us note that the phenomenological TDGL equation which models the 
overdamped (relaxational) dynamics of order parameter field to the minimum 
of the thermodynamic potential [H] is first-order in time derivative. How- 
ever, a microscopic derivation of the equation in a relativistic field theory 
using e.g., closed time path Green's function (CTPGF) leads to a second- 
order stochastic equation. Such a derivation has been done for scalar field 
theory [151 ESJ 02] • A second order TDGL equation has also been derived 
for NJL model in Ref. [18] using CTPGF method. More recently, Langevin 
equation has been derived for the chiral order parameter field in a sigma 
model including quark degrees of freedom in Ref. [19J. Here, the authors use 
influence functional method and calculate the explicit form of the damping 
coefficient as well as the form of noise correlators that appear in the Langevin 
equation which again is a second-order equation in time. This has been fur- 
ther used to discuss relaxational dynamics of the order parameter near the 
critical point [201 EI]- in condensed matter literature such a second order 
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temporal derivative term is also known as the inertial term. In this context, 
it is relevant to investigate the effect of the inertial term in the dynamical 
equation on the the late stage of ordering kinetics for the chiral transition. 
We intend to look into and investigate in which physical condition which 
evolution equation need to be considered. 

This paper is organized as follows. In Sec. [21 for the sake of complete- 
ness, we recapitulate the NJL model effective potential and the corresponding 
mapping to the M 6 Landau potential discussed in Ref. [13]. This is used to 
write down the Langevin equation for the order parameter evolution with the 
inertial term in Sec. [3j In Sec. HI we study the kinetics of chiral phase tran- 
sition with a variety of different quenches and focus on the late stage pattern 
formation dynamics. We discuss the effect of dissipation coefficient on the 
dynamics. Finally, this paper concludes with a summary and discussion in 



2 Thermodynamic potential and phase dia- 



To discuss the kinetics of chiral phase transition, we model chiral symmetry 
breaking in strong interaction in the 2-flavor NJL model. The thermody- 
namic potential in terms of constituent mass is given as [13] 
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Here, we have taken vanishing current quark mass m = 0, and introduce 
M = —2gp s , with p s = (ipip) being the scalar density and g = G[1 + 1/(4N C )]. 
We have taken the four-fermion coupling G = 5.0163 x 10 -6 MeV -2 |22j . 

Close to the phase boundary, the potential in Eq. ([1]) may be expanded 
as a Ginzburg-Landau (GL) potential in the order parameter M: 

(M) = (0) + -M 2 + -M 4 + -M 6 + 0(M 8 ) = /(M), (2) 
2 4 6 



correct upto logarithmic factors [23]. In the following, we consider the ex- 
pansion of Cl (M) upto the M 6 -term. This will prove adequate to recover the 
phase diagram in the NJL model [13J. The first two coefficients in Eq. (J2J) 
can be obtained by comparison with Eq. flTJ as 

-A 
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(3) 



We treat the higher coefficients as phenomenological parameters, which are 
obtained by fitting Q (M) in Eq. (T5j) to the integral expression for Q in Eq. (CQ). 
There are two free parameters in the microscopic theory (/i and T), so we 
consider the M 6 -GL potential with parameters b and d. For stability, we 
require d > 0. The phase diagram resulting from Eq. (j2J) is shown in Fig. [TJ 
The gap equation, /' (M) = aM + bM 3 + dM 5 = yields the three solu- 



tions for the order parameter as (M = and M\ = (—b±y/b 2 — 4ad)/2d). In 
Fig. [T]we show the phase diagram for the GL potential in [b/ (dA 2 ), a/ (dA 4 )]- 
space. For b > 0, the transition is second-order. The stationary points are 
M = (for a > 0) or M = 0, ±M+ (for a < 0). For a < 0, the preferred 

6 



equilibrium state is the one with massive quarks [23]. Next, for b < 0, the cor- 
responding solutions of the gap equations are (i) M — 0, for a > b 2 / (Ad), (ii) 
M = 0, ±M + , ±M_, for 6 2 /(4d) > a > 0, and (hi) M = 0, ±M+, for a < 0. 
As we reduce a from higher values, 5 solutions appear at a = b 2 /Ad. How- 
ever, this does not correspond to a phase transition. On further reduction 
of a, we obtain a first-order phase transition occurs at a c = 3b 2 /(16d). The 
order parameter jumps discontinuously from M = to M = ±M + , where 
M + = [3\b\/(Ad)} 1 ^ 2 . The tricritical point is located at 6 tcp = 0, a tcp = 0. In 
Fig. [U the open symbols denote 4 combinations of (b/ {dA 2 ), a/ (dA 4 )), chosen 
to represent qualitatively different shapes of the GL potential. The equation 
for the first order transition line (I) is a c = 3|6| 2 /(16<i), and for the second 
order transition line (II) is a c = 0. The dashed lines denote the spinodals 
Si and S2, with equations as 1 = and as 2 = \b\ 2 /(Ad). The typical forms of 
the Landau potential in various regions are shown in the Fig. [TJ The cross 
denotes the point where we quench the system for b < 0. 

3 The Dynamical Equation 

We intend to investigate the time-dependent evolution of the order parameter 
M(r, t) and its approach to equilibrium within the framework of the Langevin 
dynamics. The evolution of the system is described by the equation: 

^M + 7 -M + ^ = ,M). (4) 

It is based on the construction of a coarse-grained free energy functional 
Q [M] for the configuration of the order parameter field M(r, t) from a given 
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effective potential f(M) as 



Q [M] — J dr 



°L M 2 + -M 4 + -M 6 + -fvM V 
2 4 6 2 V 



(5) 



Here K measures the energy cost of spatial inhomogeneities at the domain 

interface. We substitute Eq. (J5J) in Eq. (J3J) and introduce dimensionless vari- 
ables 

M = M M f , M = ^/\aJ/\b\, 

r = £?, £ = y/Kfal 

t = tot , to = 1 / v/|o|, 

9 = \a\M 6', 7' = 7t - (6) 

Dropping primes, we then have the dimensionless form of the evolution equa- 
tion: 

— M + 7— M = -sgn (a)M - sgn (b)M 3 - AM 5 + V 2 M + 9 (f, t) , (7) 
at ot 

where sgn(x) = x/\x\ and A = |a|<i/|6| 2 > 0. The dissipation coefficient, 7 (in 
units of to) can be considered as a response coefficient that defines a time scale 
for the system. The term V 2 M arises due to spatial inhomogeneity in the 
system. The function 9 (r, t) represents a stochastic force which is assumed to 
be a Gaussian and white noise, so that (9 (r, t)) = and (9(r', t')9(r", t")\ = 
2e S(r'— r")S (f — t") [25], where e is the noise strength. In general, the actual 
form of 7 and 9 could be more complicated [26]. However, for simplicity, in 
this exploratory study we chose the simple form of Eq. (ED) which allows a 
clear distinction and comparison of the roles played by 7 and the inertial 
term in the evolution dynamics. 
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Here we will consider a system (disordered phase) which becomes ther- 
modynamically unstable when it is rapidly quenched below the critical tem- 
perature (T c ). Subsequent evolution of the system is characterized by the 
emergence and growth of domains (ordered phase). In the context of Fig. [TJ 
it corresponds to (say) a quench from a > a c to a < a c at a fixed value of 
b. The unstable homogeneous state (with M ~ 0) evolves via the emergence 
and growth of domains rich in the preferred phase (with M = ±M + ). 

In this paper, we have presented our results in dimensionless units of space 
and time. One can obtain the corresponding physical units by multiplying 
with the appropriate dimensional length-scale £ and time-scale t . For this, 
we need to estimate the strength of the interfacial energy K. We calculate 
the surface tension: o = \^K(\a\ 3 ^ 2 /\b\) J dz (dM/dz) 2 . For quark matter, a 
is poorly known and varies from 10-100 MeV/fm 2 at small temperatures [27] 
- we take a ~ 50 MeV/fm 2 . For example at T = 10 MeV and \t = 321.75 
MeV, we estimate £ = y/K/\a\ ~ 2.8 fm and t = 1/a/R - 5.1 fm [T31I28]. 

In the following, we discuss the numerical results related to the order- 
ing dynamics in the context of the phase diagram in Fig. [1] for two differ- 
ent situations. First, we study the phase-transition kinetics by quenching 
through the second-order line (II). The corresponding parameter values are 
(a/A 2 , 6, dA 2 ) = (-1.6 x 10" 2 , 9 x 10~ 2 , 7.1 x 10~ 2 ) with A = |a|d/|6| 2 = 0.14. 
This set of GL parameters correspond to (/x, T) = (231.6, 85) MeV [13]. Here 
the system evolves via spinodal decomposition. Next, we do the quenching so 
as to probe the metastable region of the phase diagram. This is achieved by 
shallow quenching through the first-order line (I) in Fig. [lj i.e., quenching to 
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the region between I and Si (marked by an arrow). Here, the system evolves 
via nucleation and growth of the droplets of favored phases. This case is stud- 
ied using the parameter values (a/A 2 , b, dk 2 ) = (3.539 x KT 3 , -0.101, 0.402) 
with A = 0.14. It corresponds to the parmeter pair (/x, T) = (322.8,10) MeV 
[T3] . This point is marked by a cross in the phase diagram of Fig. [U 

4 Kinetics of Chiral Transitions 

4.1 Case with b > 

In this case, we consider deep quenches through the second-order line (II). 
The chiral transition occurs when we quench from a > (with M — 0) to 
a < 0, where the free energy has a double- well structure (cf. Fig. [TJ. The 
chirally-symmetric phase is now unstable, and evolves to the stable massive 
phase via the spinodal decomposition. The appropriate form of the evolution 
equation is 

— M + 7— M = M - M 3 - \M 5 + V 2 M + 6(r,t) . (8) 
ot 2 at 

We solve Eq. ([8]) numerically using a simple Euler- discretization scheme with 
initial velocity dM/dt = 0. The initial state of the system i.e., M(r, 0) is 
prepared in a homogeneous state by assigning, to each lattice site, a ran- 
dom number uniformly distributed in the interval [—0.25, +0.25]. This is 
implemented on a 3-d lattice of size N 3 (N = 256), with periodic bound- 
ary conditions in all directions. For all results in this paper, we used the 
mesh sizes Ax =1.0 and At = 0.1 obtained from the linear stability analysis 
of Eq. ([8]). The thermal noise 6(f,t) is mimicked by uniformly-distributed 
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random numbers between [— A n , A n ]; where A n = 0.5 corresponding to e = 
A 2 n (Ax) d At/3 = 0.008 p3]. 

In Fig. [2] we show the evolution of an initially disordered state to Eq. OH]). 
To study the effect of dissipation, we choose 7 = 0, 0.4, and 1.0. After the 
quench, system rapidly evolves into domains of massive phase with M ~ 
— M + and M ~ M + via spinodal decomposition. The snapshots show the 
evolution of a preferred massive quark phase (M = M + ) at t — 10, 50 and 
100. The frames are the cross-section of the 3-d snapshots at z = N/2. 
As expected, the dissipation coefficient 7 controls the rapid growth achieved 
due to the inertial term in Eq. (jSJ). For 7 = 0, the inertial term is prevailing 
and therefore, we observed a rapid growth of domains of the massive phases. 
These domains are destroyed due to the oscillatory behavior of the system 
in few time steps shortly. 

However, the increasing value of 7 reduces the effect of inertial term and 
makes the system dissipative i.e., the system will now relaxed to equilibrium. 
For large dissipation (e.g. 7 = 1), the growth process is analogous to coarsen- 
ing dynamics observed in Ref . [13j • Coarsening is driven by kinks, so we can 
use the dynamics of kinks to obtain a good understanding of this evolution. 
The domain scale obeys the Allen-Calm (AC) growth law, L(t) ~ t 1 ^ 2 . In 
general, the interface velocity v ~ dL/dt ~ 1/L, where L~ l measures the 
local curvature of the interface. 

The system is characterized by a single length scale L(t) as the pattern 
morphology does not change in time apart from a scale factor. The mor- 
phology is quantitatively studied using the correlation function, C(r,t). The 
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existence of the characteristic scale results in a dynamical scaling of the cor- 
relation function: C(f,t) = g(r/L). While microscopic techniques can be 
used to measure C(f,t), scattering experiments probe its Fourier transform, 
called the structure factor, S(k,t). The structure factor is also related to 
a time-independent form or a dynamical scaling form: S(k,t) = L d f (kL), 
where g(x) and f(p) are the scaling functions jH [131 125] . 

In Fig. [31 we plot L(t) vs. t on a log- log scale. The length scale L{t) 
is obtained as the distance over which the correlation function decays to 
half its maximum value (at r = 0). The large 7 value reduces the effect of 
inertia on the system and therefore, the quenching scenario in this regime 
is similar to the one with M 4 potential (cf. Ref. [13] for the overdamped 
limit). In this case the domain growth data is consistent with the AC growth 
law, L[t) ~ t 1 / 2 at late stages. However, at smaller 7 = 0, 0.2, we observed 
an explosive domain growth with higher growth exponent. The dependence 
of the spinodal time-scale (£ sp ) over 7 is demonstrated in the inset. We 
observed a linear dependence for larger values of 7, which however, deviates 
at the smaller 7. 

4.2 Case with b < 

Let us next focus on the ordering dynamics for quenches through the first- 
order line (b < 0) in Fig. [TJ For b < 0, the first-order chiral phase transition 
occurs at a < a c = 3|6| 2 /(16<i). We consider quenches from the disordered 
state at a > a c (M = 0) to a < a c . If we quench to a < 0, the free energy has 
a double-well structure, as shown in Fig. [U Again, the ordering dynamics 
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is analogous to the previous case with b > 0. We have confirmed that the 
domain growth scenario is similar to that shown in Figs. [2H3J Subsequently, 
we focus only on quenches to < a < a c . The appropriate form of the 
dimensionless time-dependent evolution equation reads 

— M + 7— M= —M + M 3 — AM 5 + V 2 M + 9 (r, t) . (9) 
at 1 at 

The free-energy extrema of the corresponding potential: 

f(M)= l -M 2 -\M i + ±M\ (10) 

are located at M = 0, ±M + , and ±M_, where M + = [(1 + y/1 - AX) /2A] 1/2 
and M_ = [(1 - VI - 4A)/2A] 1/2 . The extrema M = 0, ±M + are the local 
minima with /(±M+) < /(0) = for A < A c . 

4.2.1 Short time behavior of the solution 

First, we study the short time behavior of the solution of Eq. (Q. We lin- 
earized it around the inflection point M by replacing M = M + <j)' and 
average over the noise. To simplify the notation we use average fluctuation 
((f>') = (f>. In Fourier space, the linearized equation becomes 

— cf>(k, t) + 7 - cf>(k, t) + [k 2 - \a\] 4>(k, t) = 0, (11) 

where a = (—1 + 3Mq — 5AMq) < in our case. Equation ( TlTl) is a homo- 
geneous second order differential equation and therefore one can write the 
general solution as 

(j)(k, t) = Ax e M ^ + A 2 e M ^\ (12) 
13 



Here A\ and A 2 are integral constants and Ai(k) and A 2 (k) are the roots of 
the quadratic equation: A 2 (k) +jA(k) + (k 2 — \a\) = 0. The wavelengths for 
which k 2 < 7 2 /4+ \a\, roots are real, that give an exponential growth in the 
system instability. Similarly, for small wavelengths, in which k 2 ^> 7 2 /4+ \a\, 
we have imaginary roots that lead to oscillatory solutions, thus preventing 
the evolution of the system. 

Further, in the strong dissipation limit (7 — > 00), the short time solution 
of Eq. dHJ» is 

<j>(k,t) = A e- (fc2 - |a|) * /7 . (13) 

Here one can see that for short wavelengths, where k 2 > \a\, the fluctuations 
4>' decay with time, while those with k 2 < \a\ grow exponentially. In this 
case the linear Eq. ( TTTj) is no longer valid and fully nonlinear Eq. ([9]) must 
be solved. However, without dissipation (7 = 0), the solutions of Eq. (fill 
become oscillatory, which prevent nucleation and growth in the system. 

4.2.2 Phase plane analysis and bubble growth 

Before we study the ordering dynamics of Eq. (Q, it is instructive to un- 
derstand the nature of the traveling-wave solutions. Our analytical under- 
standing of the domain growth problem is based on the dynamics of domain 
interfaces (kinks and antikinks) that separate regions enriched in the two 
states +M + and — M + . For the case with b > and a < 0, discussed in 
Sec. 14.11 the kinks have tanh-profiles with small corrections in the potential 
due to the M 6 -term. 
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We consider the deterministic version of Eq. (0) in d — 1 as 

— 2 M {z, t) + 1 -M (z, t) = —M + M 3 — AM 5 + —M. (14) 

We focus on traveling-wave solutions of this equation, M (z, t) = M (z — vt) 
\ I (rj) with velocity v > 0. This reduces Eq. ( fl4|) to an ordinary differential 
equation: 

(1 _ v 2) ^ + — M + M 3 — AM 5 = 0. (15) 

Eq. (JTBl is equivalent to a 2-o? dynamical system: 

dM 

y, 



dr) 



(l-v 2 ) — = M-M 3 + \M 5 --fvy. (16) 
v ' dr] 

The phase portrait of this system will enable us to identify kink solutions 
of Eq. ([H]) with these relevant fixed points (FPs) (M, y) = (0, 0), (±M_, 0), 
(±M + , 0) [13]. In Fig. HI we show the phase portraits for A = 0.14 (< A c ~ 
0.1875) and v = v s , where v s corresponds to the appearance of the saddle 
connections from — M + — > and +M + — > 0. These correspond to kinks 
traveling with velocity v s > 0, as shown in Fig. HI Here, our analysis has 
been for the case with v > 0, but it is straightforward to extend it to the 
case with v < 0. In the latter case, the portrait in Fig. H] is inverted, and the 
saddle connections (kinks) are from — > — M + and — > +M + . 

Next, we study the growth dynamics of a droplet. As we know that a first- 
order phase transition in QCD matter proceeds through a nonequilibrium 
process when the initial metastable phase breaks into droplets. In Fig. [5](a), 
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we show the growth of a droplet of the preferred massive phase (M = +M + ) 
in the background of the metastable phase (M = 0). It is obtained by solving 
Eq. Q for A = 0.14 and 9 = 0. In Fig. \5fJo), we plot v B vs. A, where v B is 
the bubble velocity. Our numerical data is in good agreement with v 8 , which 
is obtained from the phase-plane analysis. We have done a detail analysis of 
bubble growth in Ref. [T3] . 

In Fig. O the temporal variation of the order parameter field M(rj) is 
shown by solving Eq. ([16]) for v s = 0.66 at four different values of 7. In the 
absence of dissipation (7 = 0), the solution of Eq. (TIB"]) is oscillatory (solid 
line) and therefore, it is difficult for the system to escape from the metastable 
state and nucleate and grow. However, as we increase the dissipation coeffi- 
cient 7, it delays the time evolution of M(rj) to its vacuum value compared to 
the case with 7 = 0. Therefore, a large dissipation (e.g., 7 = 1) also prevents 
nucleation and growth in the system. 

4.2.3 Chiral transition kinetics 

Next, we consider the ordering dynamics of Eq. from an initial disordered 
state. We follow the same numerical scheme used in Sec. 14.11 The initial state 
with massless quarks (M = 0) is now a metastable state of the potential. 
Therefore, the thermal noise 9(f,t), which makes the system to escape from 
the metastable state have higher strength: e = 2.4. However, the asymptotic 
behavior of the domain growth in both the unstable and metastable cases is 
insensitive to the noise term [29l [30]. The chiral transition proceeds via the 
nucleation and growth of droplets of the preferred phase (M = ±M + ). 
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In Fig. we show the domain growth of the system for 7 = 0.25, 0.4 
and 0.5. The frames are the cross-section of the 3-d snapshots at z — N/2 
which show the evolution of one of the preferred phases (M = M + ) at times 
t = 20, 50 and 100. Typically, the evolution of the system begins with the 
nucleation of droplets in the early stages. The droplets grow very rapidly and 
fuse to form bicontinuous domain structures, a characteristic of the late-stage 
domain growth. The effect of 7 on the nucleation and growth can be easily 
understood by comparing the growth patterns at different 7 values in Fig. [7J 
We observed that in shallow quenching through the first order line (I), the 
system takes more time to nucleate for the limiting 7 values (i.e., 7 — > and 
7 00). However, for some moderate values (7 = 0.4), it takes less time. 

In Fig. [8]^a), we plot the scaled correlation function [C(r,t) vs. r/L] for 
the evolution shown in Fig. [TJ at three different times. For figure clarity, we 
considered only two 7 values 0.25 (open symbols) and 0.5 (filled symbols) 
respectively. At the early stages, the deviation of data sets, for both 7 
values, from the master curve reflect the morphological differences between 
the "nucleation and growth" and "domain coarsening" . As expected, scaling 
functions violate dynamical scaling in the crossover regime. However, at late 
time (t = 100) data sets collapse nicely onto a master curve for both 7 values. 
The OJK function: C (r, t) = 2/tt sin -1 exp (— r 2 /L 2 ) [31] has an excellent 
agreement with the master function. Therefore, the late-stage morphology 
in Fig. [7] is analogous to that for the unstable evolution, depicted in Fig. [2] 
(for b > 0). In Fig. [8](b), we plot the scaled structure factor [L~ 3 S(k, t) vs. 
kL] at the same times and for the same 7 values as in Fig. IH^a). At late 
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times, the tail of the structure factor shows the Porod law, S(k,t) ~ k 4 for 
k ->■ oo [32]. 

In Fig. EE we show the domain size [L(t) vs. t] on a log- log scale. For 
a comparison, we plot for different values of 7. First, we observe that the 
growth process begins once the nucleation of the droplets are over. The 
asymptotic regime is once again described by the AC growth law, L(t) ~ 
t 1 / 2 . Second, we find that the nucleation process is deeply affected by the 
dissipation coefficient 7. If we start from a moderate value of 7 the time-span 
for the nucleation increases as 7 — > or 7 — > 00. Our simulation results also 
verify the analytical results obtained by Peter Hanggi in Ref. [33], where he 
described the dissipative escape of a particle from a metastable state, known 
as Kramer's escape problem. 

5 Summary and Discussion 

To conclude this paper, we summarize and discuss the results presented here. 
We have studied the kinetics of chiral phase transitions in quantum chromo 
dynamics (QCD). At the microscopic level, these transitions are described 
by the NJL model. First, we have shown a quantitative mapping between 
the phase diagrams of the NJL model, which is an effective model of QCD 
at low energy, and the M 6 -Ginzburg Landau (GL) potential. This mapping 
enables us to identify the relevant time-scales and length-scales in QCD ki- 
netics. Near the critical point, we associated (fi,T) and the coefficients of 
the GL model. However, we have considered parameter values far from the 
critical points, and it is more appropriate to interpret the GL coefficients as 
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phenomenological quantities. 

We have studied the transitions from the massless quark (disordered) 
phase to the massive quark (ordered) phase, for different 7 values resulting 
from a sudden quench in the temperature. We have considered both cases of 
the quench - deep quench as well as shallow quench. For deep quenches, the 
massless phase is spontaneously unstable and evolves into the massive phase 
via spinodal decomposition. In purely inertial case (7 = 0), we observed a 
rapid spinidal decomposition. Domains of massive phases do not relax to any 
equilibrium value due to oscillatory behavior of the system. However, they 
do with large 7 vales. The evolution morphologies show self-similar scaling, 
and can be quantitatively characterized by the order-parameter correlation 
function or its Fourier transform, the structure factor. The domain growth 
law is L(t) ~ t 1 ' 2 . 

For shallow quenches, the massless phase is metastable. The system 
evolves via nucleation and the growth of droplets of the preferred phase. 
The dissipation factor 7 affects the time-span for nucleation and growth. We 
observed that the nucleation time does not have any significant effect in the 
late-stage of the dynamics. Once the nucleation and growth of droplets is 
over, droplets quickly merge to form bicontinuous spatial domains which is 
analogous to the late-stage domain growth; the domain growth law is the AC 
law L{t) ~ t 1 ' 2 . 

It is shown that a first-order phase transition in rapidly expanding mat- 
ter should proceed through a for-from-equilibrium process when the initial 
metastable phase nucleated into droplets of massive phases. It would be quite 
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exciting to see if QCD matters undergo such a phase transition in relativistic 
heavy-ion collision experiments by making droplets of quark-gluon plasma. 
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Figure 1: Phase diagram for the Landau free energy in Eq. ([2]) in the 
[b/(dA 2 ), a/(<iA 4 )]-plane, where we have set the three-momentum ultravio- 
let cutoff A = 653.3 MeV. A line of first-order transitions (I) meets a line of 
second-order transitions (II) at the tricritical point (tcp), which is located at 
the origin. The equation for I is a c = 3|fc| 2 /(16<i), and that for II is a c = 0. 
The dashed lines denote the spinodals Si and 5*2, with equations as 1 = and 
a>s 2 = \b\ 2 /(4:d). The typical forms of the Landau potential in various regions 
are shown in the figure. The cross denotes the point where we quench the 
system for b < 0. The point b/(dA 2 ) = 1.269, a/(dA 4 ) = —0.225 corresponds 
to the second order quench. We do not mark this point in the figure as 
it results in a loss of clarity. The open symbols denote 4 combinations of 
(b/(dA 2 ), a/(dA 4 )), chosen to represent qualitatively different shapes of the 
GL potential. 
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Figure 2: Domain evolution of a preferred massive phase (M = M + ) after 
a deep temperature quench [through the second-order line (II) in Figs. [1] 
for three different values of 7. The frames are the cross-section of the 3-d 
snapshots at z — N/2 for t = 10, 50 and 100 obtained by numerically solving 
Eq. © with a < 0, b > 0, A = 0.14. The noise strength is e = 0.008. 
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Figure 3: Time-dependence of domain size, L(t) vs. t, for the evolution in 
Fig. [2j The length scale L(t) is defined as the distance over which C(r,t) 
decays to half its maximum value. The coarsening process obeys the Allen- 
Calm (AC) growth law, L(t) ~ t 1 ! 2 for larger values of 7, such as 0.6, 0.8 and 
1.0 respectively. However, for smaller 7 = 0, 0.2 we observed an explosive 
domain growth with higher growth exponent. 
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Figure 4: Phase portraits of the dynamical system in Eq. f fT6|) for A = 0.14. 
The phase portrait for v s = 0.6690, where v s corresponds to the appearance 
of the saddle connections from —M + — > and +M + — > 0. These correspond 
to kinks traveling with velocity v s > 0. 




Figure 5: (a) Growth of droplet of the preferred phase (M = +M + ) in 
a background of the metastable phase (M = 0) for A = 0.14. We show 
the boundary of the droplet at three different times. The innermost circle 
corresponds to a droplet at time, t = 20 (b) Plot of the bubble growth 
velocity vb vs. A. The circles refer to the numerical data, while the solid line 
corresponds to the result from a phase-plane analysis. 
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Figure 6: Temporal variation of the order parameter field M(rj) as a function 
of r] for four different values of 7. Without dissipation (7 = 0), the solution 
of Eq. (Tl6|) is oscillatory [solid line] and therefore it is difficult for the system 
to nucleate and grow. 
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Figure 7: Domain evolution after a shallow temperature quench [temperature 
quench through the first-order line (I) in Figs. [1] for 7 = 0.25, 0.4 and 0.5 
respectively. The cross-section of the 3-d snapshots at z — N/2 show the 
evolution of preferred phases with M = +M + at times t = 20, 50 and 100, 
respectively. For a moderate 7, nucleation is faster (e.g., at 7 = 0.4) than 
the other two limiting 7 values. 
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Figure 8: (a) Scaling of the correlation function [C(r,t) vs. r/L] for the 
evolution shown in Fig. [7] at three different times. At early stages, the devia- 
tion of data sets from the master curve reflect the morphological differences 
between the "nucleation and growth" and "domain coarsening". However, 
at late time (t = 100) data sets collapse nicely onto a master curve for both 
7 values. The OJK function has an excellent agreement with the simulation 
result, (b) Scaling plot of the structure factor [L~ 3 S(k,t) vs. kL] for the 
same times as in (a). At late times, the tail of the structure factor shows the 
Porod law, S(k,t) ~ k~ 4 for k — > oo. 
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Figure 9: (a) Time-dependence of the domain size L(t) vs t for different 
7 values. There is no growth in the early stages when droplets are being 
nucleated. It also reflets the effect of dissipation on the nucleation and growth 
process. However, The asymptotic regime of the domain growth is consistent 
with the AC growth law, L(t) ~ t 1/2 . 
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